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ABSTRACT 

We show that the surface gravity of a transiting extrasolar planet can be calculated 
from only the spectroscopic orbit of its parent star and the analysis of its transit 
light curve. This does not require additional constraints, such as are often inferred 
from theoretical stellar models or model atmospheres. The planet's surface gravity 
can therefore be measured precisely and from only directly observable quantities. We 
outline the method and apply it to the case of the first known transiting extrasolar 
planet, HD 209458b. We find a surface gravity of g p = 9.28 ± 0.15 ms -2 , which is an 
order of magnitude more precise than the best available measurements of its mass, 
radius and density. This confirms that the planet has a much lower surface gravity 
that that predicted by published theoretical models of gas giant planets. We apply 
our method to all fourteen known transiting extrasolar planets and find a significant 
correlation between surface gravity and orbital period, which is related to the known 
correlation between mass and period. This correlation may be the underlying effect 
as surface gravity is a fundamental parameter in the evaporation of planetary atmo- 
spheres. 

Key words: stars: planetary systems — stars: individual: HD 209458 — stars: bina- 
ries: eclipsing — stars: binaries: spectroscopic — methods: data analysis 
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1 INTRODUCTION 

Since the discovery that the star HD 209458 is eclipsed 
by a planet in a short -period orbit (IHenrv et all l2000l: 
ICharb onneau et alJl200Ch it has become possible to derive 
the basic astrophysical properties of extrasolar planets and 
compare these qua ntities with theoretical predictions (e.g. 
iBaraffe et afll2003h . However, the absolute masses, radii and 
density of the transiting planet cannot be calculated directly 
from the transit light curve and the velocity variation of 
the parent star, so extra information is required in order 
to obtain them. Additional constraints can be found from 
spectral analysis of the parent star or by imposing a theo - 
retical stellar mass-radius relation (|Codv fc Sasselovll2002h . 
but this causes a dependence on theoretical stellar models or 
model atmospheres. The uncertainties in these constraints 
dominate the overal l errors in mass, radius and density (e.g. 
iKonacki et alj |2004l , limiting the accuracy with which prop- 
erties of the star and planet can be measured. 

In this work we show that the surface gravity of the 
planet can be measured directly using only the transit light 
curve and the radial velocity amplitude of the parent star. 
No additional information is required and so accurate and 
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precise surface gravity values can be obtained. As theoretical 
studies often supply predict ed values for the sur face gravities 
of planetary objects (e.g., [Baraffe et ail 120031 ) we propose 
that this quantity is very well suited for comparing observa- 
tion with theory. In addition, the surface gravity is an im- 
portant parameter in co nstructing theoretical models of the 
atmo spheres of planets |Marlev et al.| [l999: Hu bbard et al.l 
l200lft . 



After deriving an equation for surface gravity in terms 
of directly observed parameters, we illustrate this concept by 
studying HD 209458. We then apply it to the other known 
transiting extradsolar planets using results available in the 
literature. The planet HD 209458 b is known to be over- 
sized for its mass and to be strongly irradiated by the star 
it orbits. An excellent transit lig ht curve was obtained for 
HD 209458 bv lBrown etafl I^OOlf ), who used the HST/STIS 
spectrograph to obtain high-precision photometry cover- 
ing several different transit events. Pre cise radial velocity 
studies of HP 209458 are also avail able l|Henrv et al.ll2000l : 
iMazeh et aHl200d ; lNaef et alj|2004l ). 
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Table 1. Results of the modelling of the HST/STIS light curve of HD 209458. The upper part of the table gives the optimised parameters 
and the lower part gives quantities calculated from these parameters. The midpoint of the transit, T^in I > is expressed in H JD — 2 400 000. 



Limb darkening 


Linear 


Quadratic 


Square-root 


Adopted parameters 




0.12889 ±0.00042 


0.12771 ±0.00049 


0.12799 ±0.00050 




k 


0.12260 ±0.00011 


0.12097 ±0.00025 


0.12051 ±0.00037 




i (deg.) 


86.472 ±0.040 


86.665 ±0.054 


86.689 ± 0.060 


86.677 ±0.060 


To 


51659.936716 ±0.000021 


51659.936712 ±0.000021 


51659.936711 ±0.000021 




u* 


0.494 ±0.004 


0.297 ±0.027 


-0.312 ±0.127 




V* 




0.338 ±0.047 


1.356 ±0.218 




r* 


0.11481 ±0.00036 


0.11393 ±0.00042 


0.11418 ±0.00042 


0.11405 ±0.00042 


r P 


0.01408 ±0.00006 


0.01378 ± 0.00007 


0.01376 ± 0.00008 


0.01377 ±0.00008 


X red 


1.146 


1.056 


1.054 





2 SURFACE GRAVITY MEASUREMENT 



3 APPLICATION TO HD 209458 B 



The fractional radii of the star and the planet in the system 
are defined as 



R*_ 

a 



(1) 



where a is the orbital semi-major axis, and i?* and R p are 
the absolute radii of the star and planet, respectively, r* and 
r p can be directly determined from a transit light curve. 

The mass func tion of a spectroscopic binary is given by 
fe.g. lffllditchll200lh : 



f(Mp) 



(1 



2tvG 



M p J sin 6 i 
JW* + M v) 2 



(2) 



where K* is the velocity amplitude of the star, e is the orbital 
eccentricity, P is the orbital period, i is the orbital inclina- 
tion, and M* and M p are the masses of the star and planet 
respectively. Including Kepler's Third Law and solving for 
the sum of the masses of the two components gives: 



(M+ + M p ) 2 = 



2-kGMJ sin J i 



(27r) 4 a 6 



(3) 



By substituting R p — ar p into the definition of surface grav- 
ity and replacing a using Eq.|3]we find that the surface grav- 
ity of the planet, g p is given by: 



2tt (1 
~P \ 



(4) 



Eq.[4] shows that we are able to calculate the surface 
gravity of a transiting extrasolar planet from the quantities 
P, K*, e, i and r p . This can be understood intuitively be- 
cause both the radial velocity motion of the star and the 
planet's surface gravity are manifestations of the accelera- 
tion due to the gravity of the planet. A discussion in the 
context of eclipsing binaries is given by ISouthworth et all 
lj2004bl ). 

To measure g p using Eq.[4] the orbital period, P, can be 
obtained from either radial velocities or light curves of the 
system, and is typically determined very precisely compared 
to the other measurable quantities. The radial velocities also 
give e and K+, whilst the quantities i and r p can be obtained 
directly from the transit light curve. Note that it is also pos- 
sible to constrain the orbital eccentricity from observations 
of the secondary eclipse of a system. 



In order to measure the surface gravity for HD 209458 b we 
need to know r p and i. These quantities are standard param- 
eters in the analysis of transit light curves. We have chosen 
to obtain them by model ling the high-precis ion HST/STIS 
light curve presented by iBrown et al.l (|200ll ). We followed 
Brown et al. by rejecting data from the first HST orbit of 
each observed transit. 

T o model the photometric data we used the jktebop 
cod^B l|Southworth et ah 2004a ) , which is a modified versio n 
of the ebop pr ogram (|Popper fc Etzell Il98ll : lEtzell Il98ll ). 
iGimenej (|2006l ) has shown that ebop is very well suited 
to the analysis of the light curves of transiting extrasolar 
planets. Importantly for this application, jktebop has been 
extended to trea t limb darkening (LP) using several non- 
linear LP laws l|Southworth et all 120071 ). It also includes 
Monte Carlo and bootstrapping simula tion algorithms for 
error analysis (jSouthworth et al" l2004al lbh . ebop and jkte- 
bop model the tw o components of an ec l ipsing syste m using 
biaxial ellipsoids (|Nelson fc Pavislll972l ; lEtzellll975T ), so al- 
low for the deformation of the bodies from a spherical shape. 

When modelling the data we ado pted the precise orbita l 
period of 3.52474859 days given bv lKnutson et alj l|2007h . 
We fitted for the sum of the fractional radii, r*+r p , the 

T R 

ratio of the radii, k = — = -= E -, the orbital inclination, 
and the midpoint of a transit, To. We also fitted for the LP 
coefficients of the star, rather than fixing them at values 
calculated using model atmospheres, to avoid introducing a 
dependence on theoretical predictions. The linear and non- 
linear LP coefficients are denoted by tt* and respectively. 

We assumed that the planet contrib utes no light at the 
optic al wavelengths considered here (sec Wittcnmvcr et all 
2005) and that the o r bit is circular (seelLaughlin et al ] |2005l : 



Peming et al] 120051 : IWinn et al.ll2005l ). Given suggestions 



that the choice of LP law can influence the derived inclina- 
tion l|Winn et al.ll2005l ). we obta ined solutions for the lin ear, 
quadratic and square-root laws (ISouthworth et al.ll2007|) . A 
mass ratio of 0.00056 was used (|Knutson et al.l 120071 ) but 
large changes in this parameter have a negligible effect on 
the solution. 

We have calculated robust 1 a error estimates us- 



1 JKTEBOP is written in FORTRAN77 and the source code is avail- 
able at http://www.astro.keele.ac.uk/~jkt/ 
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Table 2. Surface gravity values for the known transiting extra- 
solar planets. These have been calculated using Eq.[4]with input 
parameters taken from the literature. 
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Figure 1. Best fit to the HST transit light curve of HD 209458 
using the quadratic LD law. The residuals of the fit arc shown 
offset downwards by 0.02 in flux for clarity. 



ing Monte Carlo simula tions fsee IPress et alJ Tl992l . p. 684; 
ISouthworth et all [20 04a), whic h we have previously fou nd 
to provide very reliable results l|Southworth et al1l2005al lbh. 
This procedure assumes that systematic errors are negligi- 
ble. We find no reason to suspect that significant systematic 
errors remain in the H ST light curv e after the processing of 
this data described bv lBrown et al.1 lj200l| ) (see Fig.[T]). 

The best-fitting parameters of the transit light curve are 
given in Table[T] The best-fitting model using the quadratic 
LD law is shown in Fig.[T]along with the residuals of that fit. 
The solution using linear LD can be rejected as its reduced 
X 2 is substantially larger than for the other two solutions. 
The quadratic and square-root LD laws give very similar 
solutions with reduced \ 2 values close to one. For our final 
results we adopt the mean for each parameter along with 
uncertainties from the Monte Carlo simulations (Tabled]). 
These are in good agreement with the light curve s o lution s 
obtained by ICimenezI (|2006T l and iMandel fc Ago] i|2002l ). 
both of which used the approximation that the planet is 
spherical. 

With the orbital period given bv lKnutson et ail l|2007h . 
the stellar veloci ty ampitude K+ = 85.1 ± 1.0ms -2 from 
iNaef et al. I (120041 ) . and the results of our light curve analysis 
(Table© we find the surface gravity of HD 209458 b to be 
g p = 9.28±0.15 ms~ 2 . In this case the total uncertainty in g p 
is due to almost equal contributions from the uncertainties 
in Ki, and r D . 



4 APPLICATION TO ALL KNOWN 

TRANSITING EXTRASOLAR PLANETS 

We have calculated the surface gravity values for each of the 
known transiting extrasolar planets (apart from HD 209458), 
using data taken from the literature (Table©. In several 
cases (marked with asterisks in Table© it was not possible 



* The surface gravity values for these objects have larger error 
estimates than are needed, because their fractional radii are not 
available in the literature. In these cases we have had to calculate 
them from R p and a, which are less certain than r p because of 
the need to adopt an additional constraints to calculate them 

(see text). 

References: (1) IWinn et ail (l2007d 



(2 ) iBouchv et al.l ||2005|) ; 



(3) This work ; (4) INaef et al.l j2004h: ( 5 ) IPont et al.l ll2007l); (6) 
iKonacki et"al] d2005l); (7 ) Winn et all ll2007al); (8) IPont et all 



j2004f>; (9) IWinn etaU l2007bl) ; ~fl0) lAlonso et al.l d2004l) ; 
11) Ishporer et al.l [120071 ); (12) ICameron et al.l (120071 ); (13) 
Bakos et al.l [120071) ; f 14) iMcCullough et all j2006t) ; (15) 



Sato et al.l | |2005|); (16) iGillon et al.l hOOdt); 7l7) iBouchy et al.l 
1120041); (18) iGillon et al.l (120071) ; (19) lO'Donovan et al.l j2006t) ; 



(20) ICharbonneau et al.l (120071) 

to obtain r p directly from the results available in the litera- 
ture. In these cases it had to be calculated from R p and a, 
resulting in an increased uncertainty. This is because r p is 
a parameter obtainable directly from a transit light curve, 
whereas additional constraints (for example using theoreti- 
cal stellar models) are needed to calculate a and R p . 

The orbital periods and surface gravities of all fourteen 
transiting extrasolar planets are plotted in Fig. [2] and show 
that these quantities are correlated. The linear Pearson cor- 
relation coefficient of these data is r — —0.70, indicating 
that the correlation is significant at better than the 0.5% 
level. This correlat ion is certainly related to that found by 
iMazeh et al.l l|2005l ) between the orbital periods and masses 
of the six transiting extrasolar planets then known. How- 
ever, it may be that surface gravity, rather than mass or 
radius, is the main parameter correlated with orbital period 
for these objects. Theoretical calculations have shown that 
surface gravity is a fundamental parameter in the evapora- 
tion rates of the atm ospheres of irradiated gas giant planets 
j|Lammer et al.ll2003l ). 



5 SUMMARY AND DISCUSSION 

We have shown that the surface gravity of transiting extraso- 
lar planets can be measured from analysis of the light curve 
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Figure 2. Comparison between the surface gravities and orbital 
periods of the known transiting exoplanets. Filled and open circles 
denote the systems in the upper and lower halves of Table[2] re- 
spectively. The errorbars for HD 209458 (P = 3.52 d) are smaller 
than the plotted symbol. 



and a spectroscopic orbit of the parent star. We have anal- 
ysed the HST/STIS light curve of HD 209458 (|Brown et all 
l200ll ) with the jktebop code. By combining the results 
of the light cur ve analysis with published spectroscopy 
(Nacf etal ] |2004l ) we find that the planet has a surface grav- 
ity of g p = 9.28 ± 0.15ms -2 . We stress that this measure- 
ment does not depend on theoretical stellar evolutionary 
models or model atmospheres. 

In Fig. [3] we have plotted theoretical i sochrones for ages 
of 0.5 to 10 Gyr from lBaraffe et al l (|2003h against the mass 
and surface gravity of HD 20945 8 b, adopting a mass of 
M p = 0.64±0.06 M Jup from lKnutson et all l|2007l ). The dis- 
crepancy between the observed and predicted surface gravity 
can clearly be seen. Note that these models do not include 
the effects of irradiation from the star. 

The density of a transiting extrasolar planet is of- 
ten used to compare observation with theory, but is typi- 
cally measured with a much lower precision than its sur- 
face gravity, given the same datas et. For example , the d en- 
sity of HD 209458 b derived by iKnutson etail (|2007T ) is 
345±50 kg m~ 3 . Using the mass and radius given by Knutson 
et al. leads to g p = 9.1±0.9 ms~ 2 , where the uncertainty has 
been calculated by simple error propagation. These quanti- 
ties are both much less precise and require more elaborate 
calculations than using Eq.|4] to find the surface gravity: 
g p = 9.28 ± 0.15ms -2 . 

We have applied Eq.|4] to each of the known transiting 
extrasolar planets (Tabled . The resulting surface gravities 
show a clear correlation with orbital period (Fig.[2j which is 
connected with the known correlation between orbital period 
and mass for these objects. We propose that surface gravity 
may be the underlying parameter of the correlation due to 
its influence on the evaporation rates of the atmospheres of 
short-period giant planets. 
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Figure 3. Plot of mass versus surface gravity for HD 209458 b 
compa red to the theoretical model predictions of iBaraffe et alJ 
j2003l) for ages 0.5, 1.0, 5.0 and 10.0 Gyr (from lower to higher 
logfl)' 



As g p can be very precisely measured, and can be di- 
rectly compared with theoretical models and used to con- 
struct model atmospheres of the planet, we propose that it 
is an important parameter in our understanding of short- 
period extrasolar giant planets. In the near future, the high- 
precision light curves obtained by the C0R0T and Kepler 
satellites will allow accurate surface gravity values to be ob- 
tained for the terretrial-mass transiting extrasolar planets 
which these satellites should find. 



6 ACKNOWLEDGEMENTS 

We are grateful to David Charbonneau for making the 
HST/STIS light curve of HD 209458 available on his website 
(http : / / cf a-www. harvard, edu/^dcharbon/frames .html ), 
to Pierre Maxted for discussions, and to the referee whose 
report contributed significantly to improvements in this 
work. 

JS acknowledges financial support from PPARC in the 
form of a postdoctoral research assistant position. The fol- 
lowing internet-based resources were used in research for 
this paper: the NASA Astrophysics Data System; the SIM- 
BAD database operated at CDS, Strasbourg, France; and 
the arxiv scientific paper preprint service operated by Cor- 
nell University. 



REFERENCES 

Alonso, R., et al., 2004, ApJ, 613, L153 

© 0000 RAS, MNRAS 000, 000-000 



Direct determination of surface gravities of exoplanets 5 



Bakos, G. A., et al., 2007, ApJ, 656, 552 

Baraffe, I., Chabrier, G., Barman, T. S., Allard, F., 

Hauschildt, P. H., 2003, A&A, 402, 701 
Bouchy, F., Pont, F., Santos, N. C, Melo, C, Mayor, M., 

Queloz, D., Udry, S., 2004, A&A, 421, L13 
Bouchy, F., et al., 2005, A&A, 444, L15 
Brown, T. M., Charbonneau, D., Gilliland, R. L., Noyes, 

R. W., Burrows, A., 2001, ApJ, 552, 699 
Cameron, A. C, et al., 2007, MNRAS, 375, 951 
Charbonneau, D., Brown, T. M., Latham, D. W., Mayor, 

M., 2000, ApJ, 529, L45 
Charbonneau, D., Winn, J. N., Everett, M. E., Latham, 

D. W., Holman, M. J., Esquerdo, G. A., O'Donovan, F. T., 

2007, ApJ, 658, 1322 
Cody, A. M., Sasselov, D. D., 2002, ApJ, 569, 451 
Deming, D., Seager, S., Richardson, L. J., Harrington, J., 

2005, Nature, 434, 740 
Etzel, P. B., 1975, Masters Thesis, San Diego State Uni- 
versity 

Etzel, P. B., 1981, in Carling, E. B., Kopal, Z., eds., Pho- 
tometric and Spectroscopic Binary Systems, p. Ill 

Gillon, M., Pont, F., Moutou, C, Bouchy, F., Courbin, F., 
Sohy, S., Magain, P., 2006, A&A, 459, 249 

Gillon, M., et al., 2007, A&A, in press, |astro-ph/0702i92| 

Gimenez, A., 2006, A&A, 450, 1231 

Henry, G. W., Marcy, G. W., Butler, R. P., Vogt, S. S., 

2000, ApJ, 529, L41 
Hilditch, R. W., 2001, An Introduction to Close Binary 

Stars, Cambridge University Press, Cambridge, UK 
Hubbard, W. B., Fortney, J. J., Lunine, J. I., Burrows, A., 

Sudarsky, D., Pinto, P., 2001, ApJ, 560, 413 
Knutson, H. A., Charbonneau, D., Noyes, R. W., Brown, 

T. M., Gilliland, R. L., 2007, ApJ, 655, 564 
Konacki, M., Torres, G., Sasselov, D. D., Jha, S., 2005, 

ApJ, 624, 372 
Konacki, M., et al., 2004, ApJ, 609, L37 
Lammer, H., Selsis, F., Ribas, I., Guinan, E. F., Bauer, 

S. J., Weiss, W. W, 2003, ApJ, 598, L121 
Laughlin, G., Marcy, G. W., Vogt, S. S., Fischer, D. A., 

Butler, R. P., 2005, ApJ, 629, L121 
Mandel, K., Agol, E., 2002, ApJ, 580, L171 
Marley, M. S., Gelino, C, Stephens, D., Lunine, J. I., 

Freedman, R., 1999, ApJ, 513, 879 
Mazeh, T., Zucker, S., Pont, F., 2005, MNRAS, 356, 955 
Mazeh, T., et al., 2000, ApJ, 532, L55 
McCullough, P. R., et al., 2006, ApJ, 648, 1228 
Naef, D., Mayor, M., Beuzit, J. L., Perrier, C, Queloz, D., 

Sivan, J. P., Udry, S., 2004, A&A, 414, 351 
Nelson, B., Davis, W. D., 1972, ApJ, 174, 617 
O'Donovan, F. T., et al., 2006, ApJ, 651, L61 
Pont, F., Bouchy, F., Queloz, D., Santos, N. C, Melo, C, 

Mayor, M., Udry, S., 2004, A&A, 426, L15 
Pont, F., et al., 2007, A&A, 465, 1069 
Popper, D. M., Etzel, P. B., 1981, AJ, 86, 102 
Press, W. H., Teukolsky, S. A., Vetterling, W. T., Flannery, 

B. P., 1992, Numerical recipes in FORTRAN 77. The art 

of scientific computing, Cambridge: University Press, 2nd 

ed. 

Sato, B., et al., 2005, ApJ, 633, 465 

Shporer, A., Tamuz, O., Zucker, S., Mazeh, T., 2007, MN- 
RAS, in press, astro-ph/0610556 



Southworth, J., Maxted, P. F. L., Smalley, B., 2004a, MN- 
RAS, 351, 1277 

Southworth, J., Zucker, S., Maxted, P. F. L., Smalley, B., 
2004b, MNRAS, 355, 986 

Southworth, J., Maxted, P. F. L., Smalley, B., 2005a, A&A, 
429, 645 

Southworth, J., Smalley, B., Maxted, P. F. L., Claret, A., 

Etzel, P. B., 2005b, MNRAS, 363, 529 
Southworth, J., Bruntt, H., Buzasi, D. L., 2007, A&A, in 

press 

Winn, J. N., Holman, M. J., Fuentes, C. I., 2007a, AJ, 133, 
11 

Winn, J. N., Holman, M. J., Roussanova, A., 2007b, ApJ, 
657, 1098 

Winn, J. N., et al., 2005, ApJ, 631, 1215 
Winn, J. N., et al., 2007c, AJ, 133, 1828 
Wittenmyer, R. A., et al., 2005, ApJ, 632, 1157 



© 0000 RAS, MNRAS 000, 000-000 



